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We study the three-body problem for both fermionic and bosonic cold atom gases in a parabolic 
transverse trap of lengthscale a±. For this quasi-one-dimensional (ID) problem, there is a two- 
body bound state (dimer) for any sign of the 3D scattering length a, and a confinement-induced 
scattering resonance. The fermionic three-body problem is universal and characterized by two 
atom-dimer scattering lengths, aad and bad- In the tightly bound 'dimer limit', a±/a — > cxd, we find 
bad = 0, and Oad is linked to the 3D atom-dimer scattering length. In the weakly bound 'BCS limit', 
a±/a — > — cxj, a connection to the Bethe Ansatz is established, which allows for exact results. The 
full crossover is obtained numerically. The bosonic three-body problem, however, is non-universal: 
aad and bad depend both on a±/a and on a parameter R* related to the sharpness of the resonance. 
Scattering solutions are qualitatively similar to fermionic ones. We predict the existence of a single 
confinement-induced three-body bound state (trinier) for bosons. 

PACS numbers: 03.75.Ss, 05.30.Fk, 03.65.Nk 



I. INTRODUCTION 



The physics of cold atoms has recently enjoyed a great 
amount of attention. A particularly interesting phe- 
nomenon in that context has been the experimental ob- 
servation of dimer (molecule) formation in ultracold bi- 
nary Fermi gases where a Feshbach resonance is ex- 
ploited. This has allowed experimental access to the full 
crossover from a Bose-Einstein condensate (BEG) to a 
BCS-type superfluid by simply tuning a magnetic field 
Qil^ylllilBl' Because of the Fcshbach-resonant behavior, 
the 3D scattering length a describing the s-wavc interac- 
tion strength among different fermions can be tuned al- 
most at will. For a > 0, one has a two-body bound state 
('dimer') that eventually can be Bose-condensed, while 
for a < 0, Cooper pairs are formed. The correspond- 
ing BEC-BCS crossover theory has been worked out on 
a mean- field (plus fluctuations) level 0,BS^3j Sind is 
widely believed to account for the basic experimental ob- 
servations ^3 ■ The Feshbach resonance has also interest- 
ing implications for bosonic systems, e.g., the existence of 
a similar crossover from atomic BEC to molecular BEC 
l^^^^ll^ . or Bose-enhanced 'quantum superchemistry' 

A related but different quasi-lD problem arises for ei- 
ther a two-species (t,i) Fermi gas or a single-species 
Bose gas. When such a cold atom gas is confined to 
a sufficiently tight harmonic transverse trap, it enters a 
ID regime. On the two-body level, there is always a 
bound state, even for a < 0, and one has a confinement- 
induced resonance (CIR) in the ID atom-atom scattering 
length ttaa ^^^M This 'shape' (or 'geometric') reso- 
nance is similar to a Feshbach resonance Instead 
of a magnetic field, the ratio a/a± between the 3D scat- 
tering length and the transverse confinement lengthscale 
a±, see Eq. (II. 2|) below, is now used to sweep through 
the resonance. Albeit there are strong quantum fluctu- 



ations in ID, preventing both BEC and a true BCS su- 
perfluid in the thermodynamic limit, the presence of the 
CIR leads to a rather similar scenario as for the stan- 
dard (3D) BEC-BCS crossover [H El H^. Moreover, 
this ID analogue of the 3D BEC-BCS crossover and its 
bosonic complement appear to be experimentally feasi- 
ble. Recent progress towards the realization of ID traps 
has been tremendous [Hill IHIHIlll and could lead to 
the observation of interesting aspects of ID many-body 
physics as outlined below. Moreover, on the theoreti- 
cal side, powerful many-body techniques are available in 
ID systems, e.g., bosonization or the Bethe Ansatz 
[2^|33- Such methods often allow for exact statements. 
Below we analytically solve the three-body problem for 
ultracold fermions or bosons that are confined to quasi- 
lD by a transverse trap potential. We show that, for low 
energies, atom-dimer scattering solutions can be com- 
pletely described in terms of two scattering lengths Gad 
and bad, where aad can be linked to the standard 3D 
atom-dimer scattering length. The non-standard scatter- 
ing length bad must be introduced to describe ID scat- 
tering processes in general. While the (fermionic and 
bosonic) three-body problem in 3D systems has a rather 
long history, reviewed in Rcfs. |19ll3l|. the quasi- ID sit- 
uation has not received much attention so far. 

In what follows, we assume two different hyperfine 
states of a fermion species of mass mo, or a single bo- 
son species, to be trapped in the parabolic transverse 
confinement potential 



Uc{r) = ^moivlix^ +y^), 
with associated confinement lengthscale 



a± = y^2h/moUj±. 



(1.1) 



(1.2) 



For a non-parabolic potential, the treatment is more in- 
volved |32l |. and new effects appear, e.g., additional rcso- 
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nances. Under the standard pseudopotential approxima- 
tion js^, for the low-energy limit of interest here, the 3D 
interaction can be written as 



V{r) = (5 r) — (r-), 

mo or 



(1.3) 



which assumes that the interaction range is the shortest 
lengthscale of relevance in the problem. On the two-body 
level, interactions between atoms are thereby described 
in terms of the 3D scattering length a. Normally, at low 
energy scales only s-wave interactions matter. Identi- 
cal fermions then do not interact because of the Pauli 
principle, but interactions among different ones will be 
present and are given by Eq. (|1.3|1 . The free (unconfined) 
problem is well-known to possess exactly one bound state 
('dimer') on the two-body level for a > Once the 
confinement (|1.1|1 is present, however, there is always a 
two-body bound state with dimensionless binding energy 



(1.4) 



where as = 1/«^_b is the (longitudinal) size of the dimer. 
This quantity is determined by the condition |l7l IT^ 



C(l/2,r!B)+ai/a = 0. 



(1.5) 



For our purposes, the Hurvitz zcta function can be de- 
fined via its integral representation (see appendix . 



C(i/2,r!) 



dt 



TTt V 1 - e- 



(1.6) 



Since C(l/2,r2) is monotonic in f2, there is precisely one 
bound state for any given a±/a. This bound-state wave- 
function can then be expressed in terms of the single- 
particle Green's function Ge{^, i"') for the cylindrical har- 
monic oscillator with reduced mass mo/2, see Eqs. (|2.14|l 
and l|2.15|l below. For a^^/a — oo, the 'BCS limit' is 
reached, where — {a/a±)'^ ^ 1. (For simplicity, we 
shall use the phrase 'BCS limit' also for bosons, even 
though there is no Cooper pairing in that case.) The 
dimer is then very elongated, with size as ~ 3> 
a± in the axial direction, and confined to the transverse 
ground state. The dimer in this limit can be effectively 
described by a ID contact interaction obtained by pro- 
jecting the pseudopotential H1.3|l to the lowest transverse 
state. In the tightly bound 'dimer limit', a±/a — > -|-oo, 
the dimer becomes spherically symmetric with the size 
as « a <C aj^. Here, one recovers the pseudopotential 
bound state of the unconfined problem (for a > 0) with 
the large reduced binding energy flB — (a_L/2a)^ ^ 1. 

The analogue of the Feshbach resonance in a ID con- 
fined system is then realized by the CIR. Solving the 
two-body scattering problem with just one open chan- 
nel, the ID scattering length between two atoms can be 
extracted and is found to be |0, 0| 
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C 



At low energies, this implies that one can use the ID 
atom-atom interaction potential 



2ft2 



(1.8) 



C = -C(l/2) - 1-4603. (1.7) 



The CIR, where gaa — > ±oo, then occurs for aaa ~ 0, cor- 
responding to = 1, and can be reached by tuning a±^ 
or a. The physical picture behind the emergence of the 
resonance has been elucidated in Ref. and shown to 
be similar to a Feshbach resonance. Atoms in the lowest 
transverse level (open channel) are coupled resonantly to 
a bound state in the two-body sector restricted to ex- 
cited channels (closed channel). This picture has also 
recently been extended to non-parabolic potentials |32l |. 
Although up to now, no clear experimental evidence for 
CIR behavior has been published, there are several pos- 
sibilities to observe it using standard techniques, e.g., via 
the momentum distribution or Bragg spectroscopy. 

The above considerations suggest that it is mandatory 
to analyze the consequences of this CIR on the many- 
body scattering properties [23, |^. As a first step to- 
wards a full understanding of this problem, we discuss 
the analytic solution of the corresponding three-body 
problem, which here becomes possible due to the short- 
rangedness of the atom-atom interaction. A natural ques- 
tion then concerns the scattering properties of the atom- 
dimer system, for instance, the scattering length aad- For 
fermions, this problem has been solved already a long 
time ago by Skorniakov and Ter-Martirosian (STM) |35l |. 
who found aad ~ 1.2a. (This result was recovered re- 
cently 01 as a particular limit for different masses of 
the two fermion species.) One may then wonder whether 
atom-dimer scattering will also show resonance enhance- 
ment, whether three-body physics remains universal (i.e., 
determined by two-body quantities only), and whether a 
three-body bound state ('trimer') is possible. We find 
that in the fermionic case (with equal masses), the low- 
energy physics is universal like in the 3D case |3^ . and 
there is no trimer state. However, in the bosonic case, 
important differences to these answers for fermions can 
and do arise. 

At this point, we pause in order to briefly discuss the 
bosonic three-body problem in 3D . It was shown by 
Thomas in 1935 [33 that there is no lower bound on the 
energy of a system of three bosons interacting via zero- 
range forces ('Thomas collapse'). The corresponding in- 
tegral equation for the scattering amplitude was derived 
by STM m, but in Ref. [13 it was shown that their 
equation is not well-defined, since it allows for an infi- 
nite number of solutions. A scheme of choosing the cor- 
rect solution based on the orthogonality of wavefunctions 
with different energies was suggested. A related equation 
for the bound states was investigated in Ref. where 
an infinite number of three-body bound states with ar- 
bitrarily low energies was found. Furthermore, in the 
unitary limit a ^ oo, bound states condense at the con- 
tinuum threshold. The condition of Ref. |39j was identi- 
fied as a self-adjoint extension of the involved operator. 
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Finally, Efimov solved the problem |37l | by introducing 
a three-body short-distance parameter. Away from the 
unitary limit, there is at most a finite number of bound 
states ('Efimov states'), but in the unitary limit a pe- 
culiar hierarchy of infinitely many bound states emerges 
[T^ 1^3 • This regularization also affects the scattering 
solution. To summarize, the 3D integral equations for 
bosons generally require both a long-distance cutoff (the 
scattering length a) and a short-distance cutoff. Efimov's 
real-space implementation of the short-distance cutoff is 
not very convenient for the confined problem. Fortu- 
nately, Petrov ^3 recently suggested a different and in 
the present context more useful regularization based on 
the energy dependence of the scattering length. 



(1.9) 



evaluated at the collision energy Ec of the atom-dimer 
complex, i.e., the total energy minus the (kinetic and 
confinement) energy of the relative motion. Here R* is a 
parameter related to the sharpness of the Feshbach res- 
onance, appearing also in the effective range expansion 
for the 3D scattering amplitude [33| . 



f3D{k) = - 



1 



ik + R*k^ 



(1.10) 



We shall use the regularization (|1.9|l in our discussion of 
the bosonic three-body problem. Although in strictly ID 
systems, there is no Efimov physics [l^, in the quasi-lD 
case of a confined gas under study here, Efimov states 
can and will be relevant in certain limits. 

The tunability of the scattering length a in cold atom 
experiments is obtained by using Feshbach resonances. 
For distances on the order of the interatomic poten- 
tial, the two-body problem is coupled resonantly to a 
bound molecular state involving different spin states 
(closed channel). For larger distances, there is never- 
theless only a small admixture of this closed channel and 
the two-body scattering is essentially reproduced by a 
one-channel zero-range potential with the scattering am- 
plitude (|1.1U|) . omitting the (small) background scatter- 
ing length. We refer to Ref. ^3 for more details. For 
fermions, we will assume that R* does not appear in 
Eq. (|1.1U|) under the condition that kR* <C 1. This im- 
plies in particular that R* <^ ub and i?* <C a±. 

The structure of this paper is as follows. In Sec.m we 
derive an integral equation that determines the complete 
solution of the fermionic three-body problem. This inte- 
gral equation can be solved explictly after projection to 
the transverse ground-state of the trap, as is discussed in 
Scc. lIIII The role of the higher transverse channels is then 
addressed in Sec. IIVI where also explicit contact to pre- 
vious work in the unconfined case [ssl. l36l| is established. 
A brief account of some of these results for fermions has 
been given in Ref. • The bosonic three-body problem 
is then studied in detail in Sec.0 In Sec. IVII we high- 
light the tight connections to the Bethe Ansatz existing 
in the BCS limit for fermions and bosons, and present ex- 



act results obtained from this link. Finally, we conclude 
in Sec. En 



II. FERMIONIC THREE-BODY PROBLEM 

In this section, we consider the fermionic three-body 
problem (tti) with two identical fermions, where the 
'spin' indicates two selected hyperfine states of the atom. 
We denote by Xi (x2,3) the position of the J, (the two |) 
particles, and perform an orthogonal transformation to 
variables (x, y, z) in order to decouple the centcr-of-mass 
coordinate z. 







3 









2/V3 -l/x/3 -1/V3\ /xi' 

-1 1 ^2 

'273 72/3 / \^3; 



(2.1) 



For the harmonic confinement 1)1.1(1 . the potential re- 
mains diagonal in the positions, and the Schrodingcr 
equation reads 



mo 



U,{X) - e\ ^{X) = - ^ V{r±)^{X), 
^ ± 

(2.2) 

where X = (x, y) is a six-dimensional vector. With these 
definitions, the distances between the J, particle and each 
I particle are 



r± = a/3x/2 ± y/2 = sin x ± cos Q y. 



(2.3) 



where we introduce the ang le = 7r/3 for notational con- 
venience, see also Ref. [Sg, such that •svnQ = sin(2^^) = 
73/2 and cos 6* = - cos(26') = 1/2. Equation (O then 
allows to incorporate the interactions Y[y±) via bound- 
ary conditions imposed for vanishing distances between 
t and I atoms. For r± — > 0, this implies the singular 
behavior 



*(X) ~ T 



47rr± 



(l-r±/a). 



(2.4) 



where the vectors rj^ j- = cos x =F sin Q y are orthogonal 
to r±, respectively. Since the pseudopotential H1.3|l acts 
on 'J according to 
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y(r±)«'(X) = T/(r±,±)<5(r±), 
the Schrodinger equation 1(2. 2|l becomes 



-— V|, + f/e(X)-i?U(X) = J2^^^1^S{r^) 
mo J ^ mo 

= 5(x,y). (2.5) 
Now we introduce the two-particle Green's function 

V'Ai(ri)V'A.(r2)V'Ii(r3)V'I,(i-4) 



G^'(ri,r2;r3,r4) = 



Ai ,A2 



Ex, +Ex,~E 



(2.6) 
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where ip\ denotes the eigenfunctions to the single-particle 
problem (for reduced mass mo/2) with eigenenergy Ex- 
The quantum numbers A include the longitudinal (ID) 
momentum k, the integer angular momentum m, and 
the radial quantum number n = 0, 1, 2 . . ., whence 



Ex = hu;±{2n + \m\ + 1) + ft^P/nio 
il>x = e™^i?„™(p)e*'=^ 

with radial functions [see also Eq. H1.2|l ] 
1 



Rnm{p) 



V«V(»^+ l"^|)! 



\/7r a_L 



(2.7) 



(2., 



where L™(a;) denotes standard Laguerre polynomials. 
Using Eq. H2.6|l . the general solution of Eq. (|2.5|l can be 

expressed as 

vI/(X) = vl/o(X) + j d^'dy'd^\^,y-^',y')S{^,y'), 

(2.9) 

where ^'o(X), only present for positive energy E, is a 
homogeneous (free) solution. In this paper, we restrict 
our attention to states with £" < 0, such that ^'o = 0, 
and write 



(2.10) 



where the relative momentum k of the atom-dimer com- 
plex is sent to zero later. 

Since S'(x, y) in Eq. (|2.9|l involves the variables r±,± 
and r±, see Eq. (|2.5|l . it is advantageous to switch to 
these by virtue of the orthogonal transformation 



sin 9 ± cos 9\ f X 
cos 9 =Fsin0 J \y 



(2.11) 



Similarly, we can switch from (rx__,r_) to (r^__|-, r_|-) 
using 



-cos 26* sin 261 \ / r_ 
sin 26* cos 20/ in 



(2.12) 



(2) 

Since these are orthogonal transformations, stays 
invariant. 



G^^>{^, y; x', y') = G^"' (r±, ri,±; r'±, r'^^i). 

Furthermore, the integration measure in Eq. (|2.9|) is 
also invariant, dx'dy' = dr'j^ ±d'^±- We thus find from 

Eq. ins 

9^*(r,ri) = J drl/(rl)[Gg)(r,r^;0,rl) 

- G^^ (- cos(26')r -f sin(26i)r_L , sin(26l)r + cos(26l)ri ; 0, 

where r = r_ and rj^ = Next we implement the 

r ^ limit according to Eq. (|2.4|) to obtain a closed 



equation for f{r±). This limit can be directly taken for 
the non-singular second term in Eq. (|2.13|) . while the 
first term contains the singular behavior necessary from 
Eq. (|2.4|) . Once this singular behavior is removed, one 
obtains a regular integral equation for f{r±) [s^. l3a|. 

It is then convenient to transform into the complete 
basis {V'a} specified above, f i^±) = La /AV'A(ri)- No- 
tably, the last term in Eq. H2.13|l can be expressed in 
terms of the single-particle Green's function 



Ex-E ' 



(2.14) 



since 



J drl G^')(sin(20)ri,cos(20)ri;O,rl)/(rl) = 
^G£_£,(sin(20)ri,O) i^a(cos(20) r^) fx. 

A 

The Green's function l|2.14|l for r' = has the integral 
representation 



Gi5(r,0) = 



Too 



dt 



-nt 



Vrd 1 - e- 



(2.15) 



X exp(^--^-£^coth(t/2)j, 

where = {huj± — E)/2hijj±. This can be obtained 
from the Feynman (imaginary-time) representation of the 
Green's function (|2.14() . 



/•OO 

G£(r,r')= / dtY,Mr)rx{r')e 
Jo , 



■t(Ex-E) 



and the expressions for the eigenfunctions {V'a} and 
eigencnergies Ex, see Eqs. H2.7|) and H2.8() . where only 
states with to = contribute. One has to perform a 
straightforward gaussian integral over momenta k and 
the sum over n follows from the remarkable identity 



1 



1 - e 



exp 



leading finally to Eq. (|2.15|) . Using the orthonormality 
of the {"(Aa}, the first term in the integral appearing in 
Eq. I|2.13|l can be written as 



(2.13) 



j drlG(,'^(r,rx;0,rl)/(rl) 



(2.16) 



A,Ai 



For r ^ 0, the integral representation in Eq. (|2.15|) is 
dominated by small values of t. Its asymptotic behaviour 
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is then obtained by substracting and adding the leading 
term, such that 



Too 



+ai/2M)+o{i) ), 



Gis(r,0) 



where the integral representation of the Hurvitz zeta 
function has been used. Finally, performing the 

t integral, we find that for r ^ 0, Eq. (|2.15|l has the 
asymptotic behaviour Ge{^,0) ~ {mo/4:TTh'^a±)[a±/r + 
({1/2, Q)]. This implies that the leading term in 
Eq. (|2.16|) gives TOo/(r^)/(47rfi,^r), which coincides with 
the singular part of TOo\E'(r, r±)/fi^ in the r limit. 

Using Eq. (|2.4|l . and cancelling the r ^ singular 
terms in Eq. H2.13|l . straightforward algebra leads to an 
integral equation for /(r). In the {V'a} representation, it 
reads 



(2.17) 



where we use the function, see also Eqs. (|1.6II and (|1.5|l . 

=C(l/2,f^)-C(l/2,f^B), (2.18) 
and the frequencies 

nx = nB-{a±k/2f + Ex/2hu;±. (2.19) 
We find the kernel in the form 



A 



XX' 



mo 



dr^rx{r^)iPx'icos{20)ri_) 



X GE-Ey{M'20)r^,0). 



(2.20) 



In the following, it will be more convenient to use Qb 
instead of a/a± to parametrize the interaction strength, 
see Eq. H1.5(l . 

Using Eq. H2.15|l . Ax^x' can be evaluated exphcitly, but 
before proceeding further, we shall perform a rescaling. 
Until now, / has only been considered as a function of 
the variable = [Pi^)- However, in the asymptotic 
three-body scattering solution consisting of a dimer and 
one unbound atom, the atom-dimer distance d (which is 
then much bigger than the dimer size as) coincides with 
r_L only after a proper rescaling. The asymptotic solution 
is expected to be of the form 

vI/(r,rx) = <i>o(r)x(d), (2.21) 

where r = r and the atom-dimer distance is 

d = (xi + X3)/2 - X2 = r+ - r_/2. 

Here $o(r) is the wavefunction of the confinement- 
induced two-body bound state s-^d x(d) gives the 
asymptotic solution for the scattering of the free particle 



by the dimer. The connection with / is made by looking 
at the r ^ limit of Eq. (|2.21|) . The leading term is 

vl'(r,ri)~ -^X(sin(20)rx), 
47rr 

where we have used Eq. H2.12|l to express r^. as a function 
of r and rj^. In the asymptotic limit, sin(26')r^ is thus 
the atom-dimer distance. Therefore, after the rescaling 
rj^ — > sin(26')r^, the function / matches the asymptotic 
scattering solution x- This length rescaling also implies 
wavevector rescaling, k — > k/sm.{26), as well as an ex- 
tra factor sin(20) in Ax^x' ■ In addition, from now on, 
we switch to dimensionless lengths and wavevectors by 
measuring them in units of aj_ and l/aj_, respectively. 



III. LOWEST TRANSVERSE CHANNEL 

Let us then proceed by projecting the integral equation 
(|2.17|) to the lowest transverse state (n = to = 0). The 
role of the higher transverse channels will be discussed 
in Sec. IIVI Taking into account the above rescaling, and 
noting that only to = modes have nonzero overlap with 
the lowest state, 



/■°° dk' 

mk)fk = / ^Ak,k'fk 

J-oo 27r 



(3.1) 



where 

f7fc = fis + sm'^{2e){k^ - P)/4. 
Straightforward algebra gives 



A 



k,k' 



^ ^1 + cos(46i)y 

1 



(3.2) 



p + nB + [k^ + + kk']/4: - 3fcVl6 ■ 
Note that the energy reads after the rescaling 

... ^ 'ih^ f k " 
E = -2hw±_Q.B + — 



4 Too Va_L 



(3.3) 



where k is interpreted as the relative momentum of a free 
particle with reduced mass 2too/3. 



A. Atom-dimer scattering solution 

Following STM j^, we now make an Ansatz for the 
solution of this integral equation, 



f{k) ^ 2^5{k - k) + i~f{k, k) J2 



1 



(3.4) 



with a regular function (scattering amplitude) f{k,k). 
This Ansatz gives the expected asymptotic scattering 
state after Fourier transforming to real space, 



/(z) = e*^-^ + /(sgn(2)fc,fc) e^^l^l, 



-l-oo. 
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such that standard transmission and reflection amph- 
tudes [S^l can be inferred, 

t(fc) = l + /(fc,fc), rlk) = f{-~k~k). (3.5) 

In the low-energy hmit fc, — > 0, the general expansion 
of the scattering amplitude applies, 

/(fc, k) = -1 + ikbad + ikaad + 0{k^,P, kk). (3.6) 

Here aad and had are the two ID scattering lengths for 
the atom-dimer scattering process. In general, aad and 
bad are related to the even and odd partial scattered 
waves, respectively, as will be discussed in more detail 
in Sec. IVI Al Let us just note at this point that for a 
sufficiently short-ranged potential, namely with support 
tts smaller than the typical lengthscale for the wavefunc- 
tion variations, the potential can be effectively described 
by a contact (5- interaction. This requires Us <C \aad\ and 
k <^ l/ag. In that case, odd waves are not scattered by 
the potential and hence bad = from Eq. H6.1I|I . This 
is the case in particular for the two-body problem where 
only one scattering length, Uaa in Eq. p. 7(1 . is usually 
given. 

Inserting the Ansatz (|3.4|l into Eq. (|3.1|l . we obtain 



fc2 



2fkm-k)~^vY.Jy^^m'rk) 



i \fCk,k)A,,,-k + f{-~k,~k)A,,_-, 



A 



k.k ' 



(3.7) 



where V denotes a principal value integration. One can 
check that the function C{Q,k)/{k'^ — fc^) is regular when 
k ^ k, and we used above that C{Qi.) = 0. Equation 
(|3.7(l is an inhomogeneous integral equation of the second 
kind for a given value of fc, which has a unique solution 
if the corresponding kernel is invertible. In principle, it 
may be solved numerically for any value of k to extract 
the value of aad- 

However, the subsequent analysis is simplified consid- 
erably by letting k ^ 0. Formally, we expand Eq. (|3.7|l 
in k and keep only the lowest order. To that purpose, we 
first rewrite Eq. (|3.7|l in the form 



^M^2ikf{k,k) + tV 



°° dk' - 



27rfc' 



k^-k- 

—f{k' — k, k) Af^j.,_^. 
-i[Mfc)Afe,fe + /(-fc,fc)A,,_fc 



f{k' + fc, k) ^fc,fc/+fc 



A 



k,k- 



Expanding in k and dividing by 2A;, we then get to lowest 
order: 



dk' 
27rfc' 



■jdk' Ak,k'fik',0) 



A 



k,0 



/(O, fc) + 1 
2k 



+ t; dk'Ak.k'=o- 



(3.8) 



Finally, we integrate by parts, use Eqs. (|3.2() and (|3.6|l . 
and switch to dimensionless momenta by writing k = 
2\/Q,bu. Collecting terms, we then arrive at a tractable 
integral equation for h{u) = /(m, 0). With the weakly 
r2s-dependent functions 



4-p 



G{u,u') = 

oc 



p=0 



(3.9) 
(3.10) 



this integral equation reads 
du' 



,^[G{u,u')h{u')^G{u,0)h{0)] 



/it; 



c n 



2^2 

aadV^B 



3u' 



h{u) 



(3.11) 



G{u,0) + iH{u). 



Note that the real (imaginary) part of h{u) is even (odd) 
in u. The scattering length aad finally follows from the 
real part of Eq. H3.11|l and the condition h{0) = —1, see 
Eq. H3.6|l . while bad can be extracted from the imaginary 
part of Eq. (|3.11() . The integral equation H3.11|l shows 
in particular that aad/cL± and bad/a± depend only on 
r^B, and hence only on the binding energy of the dimcr. 
This already suggests universality of the threc-fcrmion 
problem. 



B. Dimer limit 

The solution of this integral equation is discussed first 
for a±/a 3> 1, where tightly bound dimers of size as ~ 
a and large binding energy, = {a±/2a)'^ ^ 1, are 
realized. Expanding the real part of Eq. H3.11|l in l/fis, 
carefully including the ris-dependence of G{u, 0) and £, 
and using C(l/2,f] > 1) w -2V0, we obtain to first 
order 



h{u) 



3 1 



ho{u) 



41 + ^1 + 3^2/4 16 1 + 3u2/4 + + 3u2/4 

3 fis 27rii'2 \l + u^ + u'^ + uu' 1 + 
4 1 



aad 



1 



j- uu' 



UbJ \3l + u^ ^B (l + w2)s 



with G(0,0) = 4/3 and Wq = EpP(V4)^ = 4/9. Here, 
we use the zeroth-order approximation ho^u) to the full 
solution h{u). The lowest order gives 



ho{u) 



1 1 + ^l + 3u2/4 

2 l + u2 ' 



(3.12) 
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where aad = — (9/32)-\/f2B fixes ho{0) = —1. The next 
order gives from h(0) = — 1 a correction to the atom- 
dimer scattering length, such that 



aad = -Hooa^\/^ + pa^l y 
where Koo = 9/32 = 0.28125 and 



(3.13) 



9 3\/3 + 47r 3 
128 ^ 8^? 32 



0.5426. 



(3.14) 



Similarly, we can then compute h[u) to first order in 
l/Hs, and thereby obtain had from the imaginary part 
of Eq. (|3.11|l . Straightforward algebra gives 



bad/ai_ = (8/9)^!; 



'3/2 



(3.15) 



for ^ 1. The vanishing value of bad is a consequence 
of the short-rangedness of the effective atom-dimer po- 
tential ^3 ■ The support of this potential is the dimer size 
as — a and goes to zero in the dinicr limit. This validates 
a repulsive zero-range ID atom-dimer potential in the 
low-energy limit, Vadiz) = gad5{z) with gad oc {-l/aad), 
very similar to the ID atom-atom scattering potential 
(|1.8|) . In cold atom systems, the validity of our treatment 
in the dimer limit is always limited by the constraint that 
the 3D scattering length a is larger than the typical size 
of the actual atom-atom potential. 



C. Numerical solution 

Outside the dimer limit, in general a numerical solu- 
tion of Eq. (|3.11|l is necessary. We describe next how an 
accurate numerical solution to Eq. H3.11|l can be obtained 
in practice. One has to be quite careful to ensure regu- 
larity of /i(u), for which we found it beneficial to Fourier 
transform to real space, where the Fourier transformed h 
is well-behaved and allows for a quickly converging solu- 
tion of the integral equation. In order to implement an 
efficient, fast and reliable numerical solution, it is manda- 
tory to find a convenient representation of the integral 
kernel. As this is a nontrivial problem, we outline its 
solution in some detail here. 

Let us first give some auxiliary relations that will be 
useful below. The function C appearing in Eq. I|3.11|l can 
be written as 

C{u) - C(l/2,f^i3(l + 3uV4))-C(l/2,f^s) 



dt 



Jo V^l-e-*A^« 



(e-t«^*-l). 



In addition, G{u, u') in Eq. 1)3.9(1 can alternatively be 
expressed in the form 



Gob(m,w') 



dt 



1 - |e-*/OB 



^- (l+u^ +uu' +ll'^)t 



We then switch back to real space by writing 

/•oo 

h{u) = - / dze*^"g(z), (3.16) 



which leads to the Fourier transform of Eq. 1(3.1111 . The 
real part of this equation is 



dz'K{z,z')g{z') 



a_L 



'-B{z) (3.17) 



with 

K{z,z') 



Qb f du 



2 J 27rM2 ^ ' 
Ki[z,z') + K2{z~z') 



— i{z--z')u 



(3.18) 



B(z) 



^1+p/Ob|z| 



Using the above auxiliary relations, some algebra gives 
with X{t) = \z' + z/2\l^/ii and the probability function 
^{X) a convenient representation for the function Ki, 



Kiiz,z') = - 



a/3 



dt 



Jo 1 - ie-*/f2B 
Xit)<P{X{t)) + ^e-^'(*) 



(3.19) 



In a similar fashion, with Z{t) = \z — z'|/v3t, we find 



K2{z-z') 



1 /3 



dt e 



4 V TT ./o 1 - e"*/f2B 



X l^-^ + z{t)mz{t))~i] 



(3.20) 



With K{z, z') = K1+K2 and B{z) given, one numerically 
computes g{z) from Eq. 1(3. 17() and then fixes aad from 
the normalization condition J dzg{z) = 1 corresponding 
to h{0) = -1. 

Let us first discuss this program in the BCS limit, 
fig <^ 1. Importantly, the kernel K{z,z') is not in- 
vertible in this limit, since there is a zero mode solution 
go{z). This function can be found analytically in mo- 
mentum space using the Bethe Ansatz as described in 



Sec. IVI Bl with the result ho{u) = u/{l +u'^), leading 
to the real-space form goiz) = sgn(z)e~'^'. Fortunately, 
the zero mode does not affect the determination of the 
scattering length aad- To see this, note that goiz) has 
odd parity and hence does not contribute to the normal- 
ization condition J dz g{z) = 1. Since also the function 
B[z) in Eq. 1(3.17(1 is even, we can restrict the numerical 
solution to even functions g{z). On this space, K{z, z') is 
invertible even in the BCS limit, and hence the numerical 
procedure is stable and reliable. 
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FIG. 1; Scattering length aad/a,± versus dimensionless bind- 
ing energy Qb- The solid curve is the numerical solution to 
Eq. and the dotted (dashed) curves represent the an- 

alytical results in the dimer (BCS) limit, respectively. In the 
inset, bad/a.± is plotted versus Qb- The solid curve gives bad 
from the numerical solution of Eq. 13.1111 . and the dashed line 
represents Uad for comparison. 



The numerical result for aad/o.±_ as function of 51^, 
covering the full crossover, is shown in Figure Q In the 
BCS limit, where a±/a <C —1 and fls ^ 1, the real 
part of Eq. H3.11|l is solved by h{u) ~ —1/(1 + w^)^, see 
Eq. H6.20|l below, and hence we find the exact result 



aad = - a±/ 1 



0.75ai/|a|, 



(3.21) 



The numerical solution shown in Fig. ^ foi' arbitrary 
CIb nicely matches onto the analytically available limits 
H3.13() and (|3.21(l . Remarkably, around fls ~ 2.2, there is 
a zero of the atom-dimer scattering length Gad- One may 
suspect that this behavior is the analogue of the two-body 
CIR in Eq. H1.7() . However, first one may notice that the 
atom-dimer 'resonance' occurs at a different VIb than the 
CIR. In addition, it is important to check whether one 
still has a simple contact ^-interaction potential, since 
otherwise the simple relationship gad oc ^l/aad breaks 
down. 

Let us then address the numerical evaluation of bad in 
Eq. (|3.6|l . which is performed by looking at the imaginary 
part corresponding to Eq. (|3.17(l . 



dz'K{z,z')g{z') = -B2{z), 



(3.22) 



where 



■uH{u) 



du 
2^' 



p=0 



Once this integral equation is solved, bad follows as 
bad 1 



dz : 



-ig{z)] 



(3.23) 



Even though the kernel in Eq. H3.22() is the same as for 
the computation of aad, the situation is quite different in 
the BCS limit. The function B2{z) is an odd function of z 
and therefore not orthogonal to go{z). This implies that 
the integral equation (|3.22|) has no solution when ilg is 
taken directly to zero. However, for any finite fls, the 
kernel K{z., z') becomes invertible, but the corresponding 
solution of (|3.22|) diverges for fls ^ 0- The numerical 
result for bad is shown in the inset of Fig. ^ While in 
the dimer limit, bad stays small, in accordance with our 
analytical result H3.15|l . in the BCS limit, it is found to 

diverge as bad oc fig . We will see in Sec. IVII that 
this divergence originates from the reflectionless scatter- 
ing property encountered in the BCS limit. Moreover, 
these results for bad imply that one cannot use an effec- 
tive ^-potential for the atom-dimer scattering outside the 
dimer limit. Notably, the vanishing of aad for — 2.2 
does not correspond to an atom-dimer resonance since 
dad oc (— l/oad) breaks down away from the dimer limit. 



IV. ROLE OF HIGHER TRANSVERSE 
CHANNELS 

So far, the role of transverse excited states {n ^ 0) has 
been taken into account only through the exact calcula- 
tion of Ak.k' for n = m = 0. We still need to investigate 
the full Eq. (|2.17|) including higher transverse channels. 
In this section, we address their effect in detail. For the 
lowest channel n = 0, we keep using the Ansatz (|3.4|) for 
f{k). For the higher channels, we have functions /„(fc) 
instead. The full integral equation (|2.17|) [taking into ac- 
count the rescaling discussed above] leads to a system of 
coupled integral equations for f{k) and the fn{k). For 
the limit fc ^ 0, following the line of reasoning in the last 
section, we arrive at the previous integral equation for 
h{u) that now includes a coupling to the higher-channel 
modes 



hn{u) 



ifn [k = 2^yQB u \ /k 



The real and imaginary parts of {h{u),h„{u)} decouple 
in the resulting equations. For clarity, we show only the 
real part, which is sufficient to analyse the effect of the 
higher channels on the scattering length aad- After some 
algebra, we find 

— 2[G0'0(u,^.X7.') -G"'°(^.,0)/i(0)] (4.1) 

{nB[l + 3u'/4]) Hu) = ^:^^^G°'"(u,0) 

a \ 



/OO 7 / 

— G°'"'(u,^.')/in'(u'), 
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where the matrix elements G"'" {u, u') can be extracted 
from the matrix elements A^'^^, defined in Eq. H2.20|l . 
with G"'" {u,u') = UbA'^''^, and the rescaling k = 

2^/TTbu. We shall specify them in limiting cases below, 
but their general form for arbitrary parameters is of no 
interest here. Note that G^'^{u,u') = G{u,u') is given 
by Eq. (|3.9|l . Furthermore, only modes with m = to' = 
are coupled to the lowest transverse mode and thus have 
to be kept. Clearly, Eq. H3.11|l is reproduced, but now 
includes a correction due to the higher channels. The 
integral equation is then closed by 



-3/2 



du' 



[G"'^u,u')h{u') - G"'"(m,0)/i(0)] 



+ C{n + nB[l + 3uV4]) h,,{u) = 



a I r^R 



(4.2) 



The system of integral equations given by Eqs. H4.1|l and 
(14.2(1 will now be analyzed in the two limiting cases. In 
fact, we will see that in the BCS limit a±/a —f — oo, 
higher channels are completely negligible, while in the 
opposite dimer limit, they cause a renormalization of Uad 
but no qualitative change in the picture put forward in 
the last section. Moreover, in the dimer limit, we can 
solve the problem analytically and establish a connec- 
tion to the solution of the unconfined (3D) problem [s^. 
Since the efi^ect of higher channels does not cause pro- 
found changes even in the dimer limit, we conclude that 
the physical picture of Sec. lIIIl is reliable and qualitatively 
correct for all aj_/a. 



A. BCS limit 

We now show that higher channels are indeed neg- 
ligible in the BCS limit, ^ 1- As the channel 
index n has to be compared with the reduced energy 
r^B, it is intuitively clear that only small values of n 
can contribute. In addition, relevant wavevectors obey 
k, k' (X a/TTb. To make this more quantitative, we spec- 
ify A^;^,' in Eq. ifOHl for the BCS limit, which in turn 

determines G"'" {u, u') appearing in Eqs. ((4.1(1 and ((4.2(1 . 
For VIb ^ 1, the integral representation 



< n.n 

^k,k' 



^^g-{ni5+n' + [fe'^+fe2+fefe']/4}t 



X / dxe "=L„(x)L„'(x/4) 



follows from Eq. ((2.20() . Both integrals can be directly 

computed, and we find A^'^, — for n > n' within these 
approximations, otherwise there are small corrections of 



order unity. For n < n', with an rJs-independent con- 
stant C„r,', we obtain 



An.n 
^k,k' 



f^B + (fc2-h)fc'2-hfcfc')/4- 



Therefore, except for the open channel n = n' = 0, A^'^^, 
is always of order unity. Let us then analyze the scaling of 
the various terms in Eq. ((4.2() as a function of fis <C 1. 

— 1/2 

On the left hand side, the first term scales as VL^ ' , 
and the second as On the right hand side, the first 

— 1/2 

term is at most of order aad/ci± oc Jl^ , and the last 

1/2 

term scales as fig /i„. Power counting then gives hn oc 



n 



-1/2 



From Eq. ((4.1(1 . we now see that the last term 

(describing the effect of higher channels) scales as fig^, 
and is thus negligible compared to the leading terms, 
which are of order unity. This fact allows us to safely 
conclude that higher transverse levels do not affect the 
low-energy scaling behavior of aad/o.± in the BCS limit. 
This conclusion also holds for bad as one can show using 
a similar power counting reasoning. 



B. Dimer limit 

The situation is quite different in the dimer limit, 
a±/a —f +O0, where excited levels contribute to the 
asymptotic behavior of aad and bad- We first focus on 
the calculation of this correction for aad- In the dimer 
limit, relevant values for the channel number n are of 
order f2s, and we later introduce the rescaled continu- 
ous variable p = An/SflB and convert the n-summatioii 
into an integration. Keeping only leading terms in 1/D,b, 
Eqs. ((^1) and ((1?^ read 



h{u) 



aad 



4l + yT+3?IV4 a_Li 



=G"'°(u,0) (4.3) 



dp' 



^G'.p\u,urhp' in') 



and 



-2 (Vl + 3(u2+p)/4 - l) hp{u) (4.4) 



O-ad 




GP'°{U,0) 



du' 
2^ 



GP^P'iu,u')hp'{u'), 



where we use hp{u) ~ ClBhp{u). The second equation 
of this system is now independent of the first one and 
can be solved consistently. The kernel G^'^ {u, u') in the 

dimer limit VIb ^ follows from A^l'^, , sec Eq. ((2.20(1 . 
which has the integral representation 



— g-{f2f3+n' + [fc'2+feHfcfc']/4}^ 
t 



X / dxLn{x)Ln'{x/4)e "^Z*. 



(4.5) 
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To evaluate the x-intcgral, let us analyze the integral 
1^1 dx e~^^* Ln{Xx)Ln>{fix) 

DO 

dye~yLn{Xty)Ln'{lity). 



t 



Now we define A — Xnt and fl — fin't^ which stay constant 
in the limit n' — s- oo, and use an asymptotic property of 
the Laguerre polynomials, liran-too Ln{x/n) = jQ{2y/x), 
where Jq is a Bessel function. Therefore, the integral / 
in the limit of large n, n' is given by 



I ^ t dye^yja 2^ \y Ja{2^fiy 



2t 



xdxe ^ Ji 



a; ) Jo ( 2 \/Jl x 



Setting A = 1 and /i = 1/4, using the integral represen- 
tation for the Bessel function /g. 



27r 



dip 



-X cos if 



10 2^ 

and performing the i- integral in Eq. H4.5|l . we obtain 



A 



4 r^'' dip 



k.k' 



3 7o 27r 
+ 4(n + 7i' + V Tin' cos ip)/2> 



k'^ + kk')/4 
1 



After rescaling k = 2\/f^BW and n = Sp/AUb, we finally 
arrive at the expression 



GP^P {u,u') 



2tt 



dip 

2^ 



+ p + p' + v pp' cos tp 



l + u^ + u'^ + uu' 
1 



Inserting this result into Eq. H4.4|) . we can identify 
(w, </5) as the cylindrical coordinates of a 3D vector 
r (which is of course still a momentum operator), where 
h{r) is now a function of r. Writing ho{r) ~ r^/i(r), 
we obtain for the integral equation of the higher-channel 
modes (14.411 the form 



^l + 3r2/4- 1 



/io(r) + 



dr' 



hoij') 



27r2 r'2(l-Hr2 



,,12 



a-ad 



3 \ aj^ vTFj 



1 



r • r') 
(4.6) 



This integral equation is exactly the one governing the 
fermionic three-body problem without confinement |35| . 
The symmetry of Eq. H4.6|l implies that ho{v) only de- 
pends on r = |r|. We therefore write 

ho{Y) - -(2/3) (aad/a±yf^) Wo{r), (4.7) 




0.4 

0.2 ^Wi(r)\„ 




1 2 3 4 5 
r 

FIG. 2: Solutions Wo{r) to Eq. lICT . and lOi(r) to Eq. lICTI . 



where Wo{r) is the solution to 



dr' , 
In 

r r' 

Wo{r) 



l+r"^ +r''^ +rr' 



Wo{r') 



4l + ^l + 3r2/4 1 + ; 



(4.8) 



The numerical solution to this equation is shown in Fig.|21 
We find Wo(0) :^ 1.179, in accordance with STM's result 

u'o(o) w 1.2 iHiiii. 

We now go back to Eq. I|4.3|) and consider u = 0. 
Straightforward algebra then gives to leading order in 
l/r^B the atom-dimer scattering length in the dimer limit 

as 



0.636. 



(4.9) 



This value should be compared to Kqo = 9/32 = 0.28125, 
which results when higher channels are neglected, see 
Eq. (pn^ . Note that our result makes explicit contact 
to the previous solution for the unconfined case [sE HBl ■ 
This is also seen by writing Eq. (|4.9|l as 



dad = — 



red. A 



2(wo(0) a) 



(4.10) 



with the confinement scale ared.± = {3h/2mou;j_y^'^ 
for reduced mass 2mo/3 of the atom-dimer complex. 
With the 3D atom-dimer scattering length Wo{0)a ~ 
WoiO)/ Kb, this result exactly matches the dimer limit 
of the analogous two-body result H1.7|l . Equation (|4.10|) 
predicts that the product of the ID and 3D atom-dimer 
scattering lengths (in units of aj_) is a universal constant, 
independent of the atomic properties. In fact, we shall 
see later that it is even independent of statistics. 

Let us then turn to the determination of bad in the 
dimer limit, including the contribution of higher trans- 
verse channels. The calculation follows closely the one 
for aad, and we shall therefore only briefly outline the 
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various steps. Starting with the imaginary part of the 
general coupled equations resulting from H2.17|l , we arrive 
at new coupled equations very similar to Eqs. H4.1() and 
(|4.2() . We then use the asymptotic expression for A'^'^, 
in the dimcr limit and perform the rescaling k = 
n = 3p/40s and hp{u) = Qshpiu). Once again we 
identify (u, y^, (p) as the cylindrical coordinates of a 
3D vector r, and, in order to obtain an integral equa- 
tion with a spherically symmetric solution, we define 
Imft,(r) = — (it/r^) (l/r^s) wi(r). After some algebra, 
we eventually obtain an integral equation for wi (r) very 
similar to Eq. H4.8|l . namely 



reduces to 



wi{r) 



1 



4 1 + ^1 + 3r2/4 ^ Jo 
with the kernel 



dr'K{r, r')wi(r') 



1 



3(l + r2)2^ 
(4.11) 



1 



2- 



■In 



1 + + r'2 + rr' 



(4.12) 

The numerical solution to this equation is shown in Fig.|2] 
Including the higher channels, we thus find for bad instead 
of Eq. (pn5|l the result 



bad 



u.i(0) 



f7 



3/2 



u;i(0) ~ 0.952. 



(4.13) 



Notably, this is essentially the same behavior as in 

Eq. (|3.15|l . where the higher channels were neglected. 

Even the coefficient of the asymptotic behavior, 0.952 
compared to 8/9 ~ 0.889, is very close. 

To summarize, we conclude that in the deep dimer 
limit, the relevant physics is not changed by the confine- 
ment potential. 



V. BOSONIC THREE-BODY PROBLEM 

In this section, we consider the problem of three iden- 
tical bosons in a tight transverse harmonic confinement, 
allowing to analyze the ID bosonic analogue of the prob- 
lem studied up to now for fermions. The main difference 
is that the three-body wavefunction is now totally sym- 
metric. On a formal level, the calculation is similar to 
the fcrmionic one, so in what follows we go through it, 
highlighting the differences and keeping the same nota- 
tion as much as possible. The definition of the coordi- 
nate system remains the same, i.e., xi, X2 and xg are the 
3D positions of the bosons. Next the orthogonal trans- 
formation (|2.1|l to variables (x, y, z) is performed. The 
harmonic trap potential is still diagonal in the positions 
after this transformation, and the centcr-of-mass coordi- 
nate z again decouples. Hence the three-body problem 



mo 



-V2 



;7,(x)-£;U'(x 



(5.1) 



V{y) + Y.V{v±) 



*(X) 



As in Eq. (|2.3|l . y± denotes the distance between the bo- 
son at xi and the X2 or X3 particles, respectively, while 
y = X2 — X3. Note that all three bosons interact, lead- 
ing to the extra term V{y) in Eq. (|5.1|l compared to 
the fermionic Eq. H2.2|l . Incorporating atom-atom inter- 
actions via the pseudopotential approach, one boundary 
condition now reads 



*(X) 



/(x) 
47r 



while the other two are 
/(r^,±) 



*(X) 



47r 



1 



1 

a 

1 

a 



0, 



0, 



(5.2) 



(5.3) 



where r^ -t = — cos(0)x ± sin(0)y carries an extra mi- 
nus sign compared to Eq. H2.11|l . Because the three- 
body wavefunction is fully symmetric, the conditions 
(|5.3|) are redundant, and it is sufficient to satisfy only 
Eq. H5.2|l in what follows. As discussed in the Introduc- 
tion, the bosonic three-body problem requires a short- 
distance regularization R* , which we implement using 
Petrov's scheme The regularization (|1.9|l is formu- 
lated via the action of a differential operator on the scat- 
tering amplitude that describes the relative kinetic and 
confinement energy. The boundary condition (|5.2|) (for 
y ~> 0) is thereby modified to 



*(X) 



1 



R* 



moE 



/(x) 
47r '■ 



(5.4) 

where px is the component of x along the transverse di- 
rection. The boundary condition for the two-body prob- 
lem also changes, now leading to a modified equation for 
the dimensionless binding energy Q,b defined in Eq. (|1.4|1 . 
Instead of Eq. 11.5|l . it now reads 

C(l/2,f^i3) + a_L/a-4i?*f7B/aj_ = 0. (5.5) 

In the dimer limit, ils 1, the 3D result ^] follows. 



Kb 




(5.6) 



while in the BCS limit, R* gives only subleading correc- 
tions. 

In an identical way as for fermions, we then obtain the 
bosonic analogue of Eq. H2.13|l . take the y ^ limit ac- 
cording to Eq. (|5.4(l . and finally obtain the integral equa- 
tion like in Eq. H2.17|l . This modified integral equation 
for bosons reads 



/:(r!;,)/A-^AA,A'/A 



(5.7) 
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where fi^ is given by Eq. H2.19() , and the matrix elements 
are 



Ax.x = - 



dx f/'^ (x) ijjy ( — X COS 0) X 

"^0 J (5.8) 
X Gs-E'^i^ sine',0). 

When compared to the fermionic Eq. H2.20|l . the bosonic 
Eq. (|5.8(l carries an overall factor of —2, implying that 
an effective repulsion has turned into an attractive force. 
Moreover, the integrand in Eq. H5.8|l remains unaltered 
since cos 20 = — cos 9 and sin 26* = sin 6*. In addition, the 
£ function (|2.18|l is modified to 

C(nx) = C(l/2, - C(l/2, ne) - R*{k^ -P + 2n), 

(5.9) 

where from now on, all lengths (momenta) will again be 
given in units of a_L (l/aj,). Note that the distance be- 
tween the center-of-mass of the two bosons X2,3 and the 
one at xi is xsin6'. Since sind = sm(29), exactly the 
same rescaling as in the fermionic case will be employed 
in what follows. 



where the coefficient f3 is given for r* = as 

(3 ^ 9/256 + (3\/3 + 47r)/87r + 3/64 ~ 0.6638. 

We also find bad/ai_ ~ — (4/9)il^^^^, which resembles the 
fermionic equivalent (|3.15|) . 

However, as in the fermionic case, Eq. H5.10|l is not 
sufficient in the dimer limit, and higher transverse chan- 
nels must be included. Skipping details of the calculation 
- which closely parallels the fermionic one in Sec. IIVBI 
~ and keeping the same notation, we find instead of 
Eq. 14.6|l the bosonic scattering solution 



v/l + 3r2/4 - 1 + 3r*rV4 



(5.14) 



dv' 



a-ad \ 1 



a I 



/IT 



B 



1 



Defining the function Woir) as in Eq. (|4.7|) . we obtain 
after angle integration instead of Eq. H4.8|l the equation 



A. Scattering solution 

Let us then proceed by projecting Eq. (|5.7() onto the 
ground state. We use the same scattering Ansatz (|3.4|l 
as for fermions, and after expanding in fc = 2\J^b u, one 
easily obtains the bosonic version of Eq. (|3.11|) , 



du' 

[G{u,u')h{u')-G{u,0)hm 

2 27rit'^ 

^""CiuMu) = ^M^/^G(u,0)+*i7(u) 



(5.10) 



4u2 



a± 



with the functions G{u,u') and H{u) defined as for 
fermions, see Eqs. (|3.9|) and (|3.10|) . Furthermore, 
Eq. H5.9|l gives 

Ciu) - C(l/2,f^s(l + 3wV4)) (5.11) 
- Cil/2,nB) -3R*nBU^. 

In the bosonic case, the atom-dimer scattering length 
o-ad/o-L therefore depends on the two dimensionless pa- 
rameters VLb and R* /a \ . In that sense, the bosonic prob- 
lem is non-universal [13 . 

Let us start with the dimcr limit, ils ^ 1. In that case 
it is useful to introduce the dimensionless regularization 
parameter 



KbR* 




(5.12) 



The solution to Eq. H5.10() in the dimer limit can be found 
analytically again, with the result 



a I 



9 



64(1 + 2r* 



nil 



(5.13) 



v/l + 3r2/4 - 1 + 3r*rV4 



dr' 



2r' 



In 



Wo{r) 



l + r^ +r'^ + rr' 
1 + + r'2 — r r' 



Wo(/).(5.15) 



At r* = 0, this is precisely STM's Eq. (15) [33, once 
the 3D scattering length is identified with Wo{0)aB = 
Wo(0)/ Kb- For finite r*, it is equivalent to the equation 
recently studied by Petrov, since Eq. (|5.15() follows from 
Eq. (12) of Ref. |4l| upon substitution of the standard 
3D scattering Ansatz. The scattering length can be ex- 
tracted in a way similar to the fermionic case. For that 
purpose, we put u = in the bosonic version of Eq. (|4.3|l . 
and arrive at 



--(1 + 2r*) = 

8 3 a^^/flB 



TT Jo 1+^2 



On the other hand, taking r ^ in Eq. H5.15|l yields 



^(1 + 2.*)^ 



Consequently we have 
aad 



2 

-1 + - 



°° , Mr) 
dr ;r 



4u;o(0)' 



(5.16) 



i.e., the same relation as for fermions, leading to 
Eq. H4.10|l and matching the analogous two-body result 
(|1.7|) . We conclude that in the dimer limit, for both 
fermions and bosons, the ID scattering length Gad is al- 
ways inversely proportional to the 3D scattering length 
Wo{0)Jkb- The latter has been studied in detail in 
Ref. kll as a function of R* /a. It was found to diverge 
whenever a new Efimov state splits from the continuum, 
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going through zero in between the resonances, see Fig. 1 
of Ref. According to Eq. (|5.1()|) . the atom-dimer 

scattering length aad will behave in the opposite manner, 
i.e., it vanishes every time a new bound state emerges 
and diverges in between. Note that in order to derive 
Eq. (|5.14|l . one needs to impose Wo{0)/kb < a±, which 
follows from ^ 1, so that divergences of Wo{0)/kb 
will in practice be smeared out on lengthscales « a_L 2> a. 

In the limit of large R* , one can solve Eq. 115.1511 ana- 
lytically by expanding in inverse powers of r* [4l|, 



Wo{r) 



3r* 1 



0{l/r 



*2\ 



Taking into account Eqs. (|5.6|l and H5.12|l . it is easily 
seen that the 3D scattering length becomes independent 
of R* , while for the ID scattering length, the contribution 
of higher channels clearly becomes negligible. We then 
obtain from Eq. (|5.13|l 



- 

" 64 a 



(5.17) 



which is a universal result in the sense that it docs not 
depend on R* . 

Higher channels can also be taken into account exactly 
for the calculation of bad in the dimer limit. We follow 
the fermionic case in Sec. II V Bl and obtain 



bad 



3/2 ' 



(5.18) 



where, in contrast to the fermionic case, wi{r) solves the 
integral equation 



Wi{r) 



'•*Wi{r) 



4l+ ^1 + 3r2/4 8 

f°° dr' 1 

Jo TT 3(1 +H)^ 



(5.19) 



with K{r,r') given by Eq. (|4.12(l . The solution for wi(0) 
as a function of R* /a is shown in Fig. |31 Apparently, 
there is no divergence due to Efimov states. In particular, 
the R* ^ limit is well-defined, with 'Wi{0) ~ 1.59, and 



-12.7aVai> 



(5.20) 



validating a zero-range atom-dimer potential in the low- 
energy limit. The other limit is R* 3> a, where wi(0) ~ 
(8/9)-\/a/i?*, and therefore 

bad = -(64/9)i?*(a/ax)'. 



(5.21) 



Comparing this result to Eq. H5.17|l . we find that the zero- 
range potential approximation requires R*a? /a\ <C 1, 
which self-consistently holds for VLb ^ 1. 

In the BCS limit, Vlg <C 1, the same argument as in 
Sec. IIV Al can be applied to the bosonic case. Therefore 
higher channels can again be disregarded in that limit. 




0.1 1 
R*/a 



FIG. 3: Solution mi(0) (full line) to Eq. ifCT^ as a function 
of R* /a. The dotted curve gives the asymptotic behavior 
tt;i(0) = (8/9) ^a/i?* for large R' /a. 



Moreover, the regularization parameter R* drops out to 
leading order, since relevant momenta are small (oc \/TIb) 
and the short-ranged cutoff ceases to matter. In this 
sense, the three-body problem becomes universal again 
in the BCS limit. The asymptotic behavior of aad and 
bad can be analytically computed using the Bethe Ansatz, 
see Sec. IVII One can verify that h{u) ~ — 2u/(l 4- v?Y 
solves the imaginary part of Eq. (|5.1U|I in the BCS limit, 
leading to 



bad 



h'{0) 



1 



(5.22) 



For the real part of Eq. (|5.1U|) , we obtain from the Bethe 
Ansatz the solution 



h{u) = 4- 



1 



(5.23) 



in the BCS limit. Evidently, the condition /i(0) = — 1 
cannot be fulfilled for any finite aad- Similarly to bad 
for the fermionic case, the asymptotic behavior of aad 
is therefore expected to diverge as aad oc — 1/fi'^''^ when 
approaching the BCS limit. We shall discuss this point 
in detail in Sec. IVII 

For the general case of arbitrary aj_/a, we have solved 
numerically the real-space version of Eq. H5.10|l as in 
Sec. nil CI neglecting higher channels. For the sake of 
simplicity, we take R* = 0. When compared to the 
fermionic Eq. H3.17|l . the kernel l|3.18|l is now K = 
Ki - A'2/2, with Ki^2 given in Eqs. and ifT^ . 

respectively. The results are shown in Fig. 2J Clearly, 
they closely match the predicted behavior in both limits. 
A finite R* in Eq. (|3.17|) gives similar results. From the 
above discussion, this is clear in the BCS limit. More- 
over, when just keeping the lowest transverse channel, 
the qualitative behavior is not modified in the dimer limit 
either. Therefore we expect a very similar picture for ar- 
bitrary R* as the one shown for i?* = in Fig. ^ as 
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FIG. 4: Plot of ttad and bad (full lines) as function of the 
dimensionless binding energy Hb for the three-boson problem 
at R* = 0. The asymptotic behavior is shown (i) for bad in 
the BCS limit, — I/^/TTb, and (ii) for aad in the dimer limit, 
(9/64) VTTb + 0.6638/ V^is. 



long as higher channels can be neglected [as presumed 
in Eq. (|3.17() ]. While this approximation is justified in 
the BCS limit, in the dimer limit it does not capture 
the correct asymptotic behavior, except for large R* /a, 
see above. In particular, the influence of Efimov states 
on aad becomes important and affects the physics in the 
dimer limit in a significant way. This is in contrast to 
the fermionic case, where the projection onto the lowest 
transverse state already yields the qualitatively correct 
behavior even in the dimer limit. 



B. Bound States: Trimers 

For the bosonic problem, one expects atom-dimer 
bound states (trimers) to form under certain conditions. 
We therefore next derive the relevant integral equa- 
tion. In this case, the scattering Ansatz cannot be used, 
and the resulting homogeneous integral equation must 
be considered. The total energy of the system is now 
E = —2 hwiSls — Et, where the trimer binding energy 
is written as 



3ff 
4 mo' 



-Kr>U. 



(5.24) 



with uq being the inverse size of the trimer in units of 
as = Projecting Eq. H5.7|l onto the transverse 

ground state, we find 




f{u) 

du'Guo "')/(w')) 



(5.25) 



where the function G,ij,(u,w') is defined by 

1 -I- 37/2 /4 _^ ,,2 

p=0 



G 



3u§/4 + v? + u'2 + uu' + p/VLi 



Equation (|5.25|) is then an eigenvalue equation for mq. 

First, we study Eq. H5.25|l in the dimer limit. In order 
to obtain a nontrivial solution, we need to consider ito <C 
1. Then f{u) is dominated by its small-it behavior, and 
using Eq. (|5.12() , Eq. (|5.25() is simplified to give 



f{u) 



32 



^2 97rOi3(l + 2r*) 



du'f{u'), (5.26) 



which implies f{u) oc l/{u'^^u'^). Self-consistency of this 
expression leads to 



32 



9f7B(l + 2r*) 



(5.27) 



Note that we have supposed here that f{u) is an even 
function. Taking instead an odd function, we get 



32 



1 



9 Trnsil + ^r*) 



du'u'fiu'), 



which does not allow for a non-trivial and self-consistent 
solution. Therefore f{u) cannot be odd in the dimer 
limit. For large r* , the trimer binding energy becomes 
universal, and is given by 



1024 fi^ 
27 Too a\ 



(5.28) 



Note that the existence of this novel confinement-induced 
trimer state (CIT) is consistent with the positive ID 
atom-dimer scattering length for bosons found in the 
dimer limit, see Eq. (|5.17|) . 

In the BCS limit, Eq. H5.25|l for the bound states re- 
duces to 



1 



1 



2 du'f[u') 

TT 



1 -f iul/A + u'^ + u'2 uu' 



(5.29) 



The fact that one can use the Bethe Ansatz to solve di- 
rectly the BCS limit also holds for bound states. Indeed, 
we will verify in Sec. IVII that Eq. 1)5. 29|) follows from a 
ID Schrodinger equation, see Eq. (|6.4|l below, and leads 
to the exact solution 



/(") 



1 



(5.30) 



for uq — 2. We have verified numerically that there are 
no other solutions, neither even nor odd. In fact, in real 
space, Eq. (|5.3U|I is nothing but the known bound-state 
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FIG. 5: Parameter uo appearing in the trimer binding energy, 
Eq. 15.24L as a function of f2s for various values of R* /a±. 
The dotted lines correspond to the asymptotic behavior in the 
dimer limit, Eq. (15.271 . 



FIG. 6: Reduced energy e = (ET/huji_){aj_/af for the GIT 
in the dimer limit as a function of R* /a. 



Bethe function for three attractively interacting bosons 

Ip CX exp [-(|zi - Z2I + \Z2 - Z3I + |Z3 - Zi\)/aaa] ■ 

(5.31) 

In between the BCS and the dimer limits, we have 
investigated Eq. H5.25|l numerically, see Fig. [S] Quali- 
tatively, we see that the trimer energy barely depends 
on R* , and in the BCS limit even becomes completely 
independent of this regularization parameter. However, 
Eq. H5.25|l is only approximate since it is derived by pro- 
jecting Eq. 1)5. 7|l onto the lowest transverse state, and 
therefore does not include the effect of higher channels. 
Notably, Efimov trimer states cannot be recovered within 
this approximation. Such states are irrelevant in the BCS 
limit, where the lengthscale Uaa of the confined trimer 
state [see Eq. (|5.31|) ] is larger than the typical length- 
scale a of Efimov states, so that the two problems decou- 
ple. However, when moving towards the dimer limit, the 
interplay between the confined trimer state and Efimov 
states plays a role. 

In the dimer limit, it is possible to study the CIT in a 
quantitative manner, including the way it is affected by 
Efimov trimer states. In fact, we have shown in Sec. lV Al 
for the dimer limit, that the low-energy atom-dimer scat- 
tering properties are described by a ID contact potential 
V{z) = (— 2fi,^/moacirf)(5(z). This potential has a bound 
state for aad > and its energy is given by 

Et = — (5.32) 

The corresponding ID scattering length Uad is related 
through Eq. (|4.1U|) to the 3D atom-dimer scattering 
length Wo{0)aB (i-G-, without external confinement) ob- 
tained from the solution of Eg. I| 5.15() . This quantity is 
shown, e.g., in Fig. 1 of Ref. [IJ. The bound state of 
this ID contact potential - present only for aad > ~- 
describes the CIT as long as its energy Et < huj±. This 



implies aad > which in turn is equivalent to the con- 
dition Wo{0)aB < a±. This coincides with the condition 
in Sec. IV Al for the validity of Eq. (|5.14() . and therefore 
of Eq. H4.10|l for bosons. The CIT is thus present only 
for Wo{0)aB < 0, and its energy follows from Eqs. (|5.32() 
and H4.10|l in the form 

^ 32 /MOWy^ (5.33) 

nuji^ \a± J 9 \ a ) 

From the numerical solution of Eq. (|5.15|l , we obtain re- 
sults for Et as a function of R* /a^ see Fig. El Obviously, 
in the dimer limit, the CIT only exists for certain val- 
ues of R* /a. While it is present for all values larger 
than R* /a ~ 2.2, below this threshold, it disappears and 
reappears periodically for decreasing R*/a. One such 
cycle is shown in Fig. El where the CIT has disappeared 
in the window 0.17 < R*/a < 2.2, but reappears for 
R*/a < 0.17. For R* < a, the behaviour is exactly 
periodic as a function of ln(i?*/a) and can be obtained 
analytically using the known result |l9j | 

= 1.46 -2.15 tan [1.006241n(aA*) + 0.09], 

(5.34) 

together with A* ~ 6.6/ R* jij. The ultimate fate of the 
CIT and its evolution in between the BCS and the dimer 
limits remains a difficult open question that is outside 
the scope of our paper. 

VI. BETHE ANSATZ IN THE BCS LIMIT 

In the BCS limit, Hb ^ 0, the role of transverse 
excited states is negligible. Confined in the transverse 
ground state, the motion of particles becomes effectively 
ID, with contact ^-interactions between bosons or be- 
tween distinguishable fermions. This strong simplifica- 
tion of the original 3D hamiltonian hence allows for ex- 
act solutions on the many-body level, namely the Lieb- 
Linigcr model for bosons [4J|, or the Yang-Gaudin model 



16 



for fermions |45j, both based on the powerful Bethc 
Ansatz. Recently, these two models have been used to 
study many-body properties of cold atoms in ID confined 
geometries |23, Hj. In this section, we use the Bethe 
Ansatz to solve the three-body problem for bosons 
and fermions in the BCS limit. This leads to analytical 
predictions for aad and bad, as well as for the bosonic 
trimcr state. We thereby derive some of the results pre- 
sented in previous sections. Noting that the atom-dimer 
scattering does not lead to any reflected wave, we also 
clarify the connection between the divergence of Uad (bad) 
for bosons (fermions) and the fact that the scattering 
process is refiectionless. The reduction of the original 
problem to a ID model is not completely straightfor- 
ward and needs to be considered cautiously. Even in the 
BCS limit, the few-body wavefunction is not restricted 
to the transverse ground state. In order to build up a 
boundary condition like Eq. (|2.4|l . imposed by the two- 
body potential for small particle distances, many excited 
transverse states have to be involved. However, as soon 
as distances between particles are larger than a±, essen- 
tially only the transverse ground state can be occupied. 
In our case, the relevant lengthscale is the ID scattering 



length Oq 



/2\a\ a_L, which justifies the neglect 



of transverse excited states. This simplification becomes 
incorrect for problems involving lengthscales < a, e.g., 
when describing deeply bound Efimov states or three- 
body recombination processes. 

Let us first explicit the connection between our integral 
equations in the BCS limit and the ID reduced problem. 
We start with bosons. The ID Schrodinger equation is 



2moE 



Ipizi, Z2, Z3) 



[S{zi - Z2) + 5{zi - Z3) + S{Z2 - Z3)] ■4!{zi,Z2, Z3), 

(6.1) 

where the two-body ID contact interaction can be recov- 
ered by projecting the 3D pseudopotential to the lowest 
transverse state, 



-S{z) = [ pdpd(t> Rooip)^^^^6{r)-^{r-)Roo{p). 
J mo or 

(6.2) 

We then closely follow Sec. |^ with the ID boundary 
condition for y 0, 



^(x,y) = /(x) 1-^ , 



(6.3) 



obtained by projecting Eq. H2.4|) to the lowest transverse 
state. One finally obtains, after some algebra and the 
rcscaling k = {\/3/2 aaa)u, the integral equation 



^l-f 3(u2 - w2)/4 
2 



/(«) 



du'f{u') 



(6.4) 



which can also be obtained from the integral equation 
(|5.7|) in the BCS limit, i.e., from the scattering approach. 
For fermions, the same line of reasoning can be followed, 
and one obtains the integral equation 16.4|l with an over- 
all factor —1/2 factor on the right hand side. Wc con- 
clude that the direct solution of Eq. (|6.1() can provide 
exact results in the BCS limit for aad, bad, and possible 
trimer bound states. 



A. Bosonic case 

Equation H6.1|l can be solved by using the Bethe 
Ansatz. The corresponding bosonic three-body wave- 
function is written in the fundamental domain Vi = 
{zi < Z2 < Z3} as 

3 

i^{zi,Z2, Z3) ^NY,A{P) e'^-i'''"' , (6.5) 
p 

where is a normalization constant, and the sum ex- 
tends over all permutations of {1,2,3}, with coefficients 
l2i 



A{P) 



n 

1<J<'<3 



2i/aa 



kpi - kpj 



(6.6) 



In other domains, the wavefunction is recovered by sym- 
metry arguments. There is only one trimer bound-state 
solution EElil, 



Ipizi, Z2,Z3) = A/'exp [-2(Z3 - Zi)/aaa] , 



(6.7) 



here given in the fundamental domain or in general 
form by Eq. (|5.31() . Using the boundary condition H6.3() . 
we find that f{u) 1/(^2 + 4) solves Eq. H6.4|) with 
u — iuQ ~ 2i, as discussed in Sec. IV Bl 

We now turn to the atom-dimer scattering problem, 
putting Uaah = u, Uaah = -i-u/2 and aaakz = i-u/2, 
with the corresponding total energy 



E ■ 



mo ala 



-1 



(6.8) 



1 + u2 -f it'2 + uu' - 3m2/4' 



see Eq. (jnSll in the BCS limit. From Eq. ^i^, the three- 
boson wavefunction is explicitely given. Using then the 
ID boundary condition 16.3|l . we find in real space after 
some algebra, up to an overall normalization constant, 

fix) = {2^3ieu) [-(2 - ie?l)e""^ + 2(2 + k?l)e-™^/2g-|x| 

(6.9) " 

where e = sgn(a:), and we have performed the rescaling 
V3x/2aaa ^ a; in order to give the correct atom-dimer 
distance in units of aaa- Remarkably, the atom-dimer 
scattering does not give any reflected wave. The atom 
just passes the dimer and only acquires a phase shift, 
without being actually backscattered. Consequently, one 
can neither use the scattering Ansatz (|3.4|l nor define Uad 
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and bad as in Eq. H3.6|l . We therefore need a more general 
definition to take into account the reflectionless situation 
found in the BCS Hmit. From Eq. H6.9() and its counter- 
part with u — > —u, we can use a simple but convenient 
basis for the incoming waves. The symmetric (antisym- 
metric) state is taken as a plane wave e™^ coming from 
the left plus a plane wave e~™^ ( — e~*"^) from the right. 
In that basis, the scattering process is described by the 
2x2 scattering matrix [sof 



S 











t + r 
t-rT 



(6.10) 



where r and t are again energy-dependent reflection 
and transmission coefficients that can be read off from 
Eq. Ht).9|) . At low energy, after expanding in fc = u/aaa, 
the phases da,s define two scattering lengths according to 

<5,(fc) = -2kaad + 0{P), SaCk) = 2~kbad + 0{W). 

(6.11) 

This definition is more general than Eq. and re- 

duces to it when the Ansatz (|3.4|) applies. We also ob- 
serve that in general two scattering lengths are needed 
to completely describe an arbitrary scattering process in 
ID. 

Using Eq. (|6.9|l . we find for this atom-dimcr scattering 
problem the result bad = — 2aaa = — l/V^s- This agrees 
with the asymptotic result in Fig. 0] for the BCS limit. 
Furthermore, we find aad ~ — oo, in accordance with the 
divergence aad oc -l/^^^''^ see Fig. H Finally, we infer 
analytical results from Eq. (|6.9() for the integral equation 
(I5.1U|) in the BCS limit. Switching to Fourier space and 
expanding / in iZ as f{u) = fo{u) + ufi{u), we find 



fa = -716(11) + 
fi = -7tS'{u) + 



1 + ^2' 

2 



m(1 4- u^) 



2\2 ' 



(6.12) 
(6.13) 



where each of these functions satisfies Eq. (|6.4|l for u — 0. 
Plugging fi into Eq. (|6.4() for u = 0, we confirm the 
asymptotic behavior bad = — l/V^B in the BCS limit, 
see Eq. (|5.22ll . Plugging /q into Eq. (|6.4|) for S = 0, we 
obtain the solution H5.23|l with h{Q) = for any finite Uad- 
Therefore the condition h{0) = — 1 implies aad ~ —00. 



B. Fermionic case 

We now consider the fermionic (tti) three-body prob- 
lem and solve it via the Bethc Ansatz. The Schrodinger 
equation is essentially Eq. (|6.1|1 with the three (5-bmctions 
replaced by (5(zi — Z2)-l-(5(2;i — Z3). This equation does not 
lead to any trimer state. The Bethe Ansatz for fermions 
is slightly different from the one for bosons. The form 
Eq. H6.5|l also applies to the fermionic case, but the coef- 
ficients A(P) are not given by Eq. H6.6|) anymore. More- 
over, the fundamental domain Pi is not sufficient, and 



we need to know the wavefunction in one more domain 
[isj . e.g., 2?2 = {z2 < zi < Z3}. The wavefunction in 
other domains then follows by anti-symmetry properties. 

To study atom-dimer scattering, we now consider the 
momenta aaaki = —i — u/2, aaak2 = i — u/2, and aaaks = 
M, with total energy H6.8|l . Imposing the fermionic anti- 
symmetry and normalizability, the problem reduces to 
the determination of three variables q;i,2,3- The wave- 
function is given in Vi by 

(6.14) 

and in V2 by 



i(k2Zl+kiZ2+k3Z3) 



(6.15) 



The variables are then obtained by imposing boundary 
conditions for the wavefunction and its first derivative 
with respect to zi and Z2 at the boundary between I?i 
and I?2, {zi — Z2 < 23}. The result is ai = 1 — x, 0^2 = x 
and as = —1, where 



X 



3m 



Once the wavefunction ^/;(zi, Z2, 23) is known, we can 
use the boundary condition 16.3|l and the rcscaling 
^/3.T/2aaa X to get 



/(x) ~ e(2 — 3iew 



g iux/2^-\x\ 



(6.16) 



where again e = sgn(a:;). Similar to the bosonic case, for 
arbitrary w, there is no reflected wave. The scattering 
only leads to a phase shift that goes to tt at zero energy. 
This implies a transmission coefficient i = — 1, in contrast 
to the bosonic result t — 1. This discrepancy originates 
from the antisymmetry of the fermionic wavefunction 
with respect to exchange of the two f atoms. As has been 
discussed in Sec. lIII CI we obtain from Eq. (|6.16|) and the 
definition H6.11|l that aad — {il'2)aaa = 0.75/VITb, and 
bad ~ +00. These results confirm the asymptotic behav- 
iors shown in Fig. ^ for the BCS limit. 

Also in the fermionic case, we can obtain analytical 
results for the integral equation (|3.11|) in the BCS limit. 
We first Fourier transform Eq. 16.16|l . The lowest order 



m u gives 



/o(") = 



u{l +u2)' 



(6.17) 



Plugged into the fermionic equivalent of l|6.4|) for u ~ 0, 
one finds the zero mode responsible for the divergence 
of bad, see Sec. IIII CI More precisely, ho = {—i)v? fa{u) 
satisfies the imaginary part of Eq. H3.11|l in the BCS limit. 
The next order in u cannot be taken directly, since the 
resulting integral in the fermionic equivalent of Eq. H6.4|l 
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0. Instead, we use 
1 1 



is not well-defined for w and u 
that f{u) has the following form 



f(u) = /o(w) + — iiTuSiu — u) + iT(u, u) 

2 \u — u u ^ 

(6.18) 

where T{u, u) is a regular function. We insert this expres- 
sion into the fermionic equivalent of Eq. (|6.4|l and look 
at the u ^ limit. After some algebra as in Sec. 1111 Al 
we obtain 



du' 
27ru'2 

1 



[G{u,u')h{u') ^ G{u,0)h{0)] 



^1 + 3u2/4 



with 

hiu) =r(M,0) 



h{u) 



(1+1*2)2 



(6.19) 
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(1+^2)2 



(6.20) 



and G{u,u') = (1 + ?i2 + u'^ + uu')^^ . This shows ex- 
plicitly that Eq. (|6.2U|) is a solution to the real part of 
Eq. l|3.11|l in the BCS limit, with the expected asymp- 
totic behavior aad/o,!. = 0. 75/ VITb. Although the form 
(|6.18() seems similar to the general Ansatz (|3.4|) , it is dif- 
ferent since there is no reflection in the BCS limit, and 
Eq. H3.4|l cannot apply. 



C. Reflectionless potential 

We have seen above that the divergence of bad (oad) 
for fermions (bosons) in the BCS limit was due to the 
reflectionless character of the atom-dimer effective po- 
tential. We shall now generalize this idea and show that 
for ID scattering processes of a particle by any (possibly 
non-local) symmetric potential, the following statements 
are equivalent: (i) there is a divergence of the ID scat- 
tering length ai (the 'odd' scattering length bi), (ii) the 
potential has a quasi-bound state with even (odd) parity 
precisely at zero energy, and (iii) there is no reflection 
at zero energy. The atom-dimer scattering problem in 
the BCS limit, where higher channels are negligible, is in 
fact equivalent to a ID scattering problem with such an 
effective non-local potential ^^l- Therefore, the consid- 
erations in this section directly lead to physical insights 
for the BCS limit. The fact that in our case the reflection 
coefficient is zero for any incoming energy is an additional 
feature that cannot be inferred from the behavior of ai 
and/or bi. 

Symmet ric p otentials imply the low-energy scattering 
properties 



ricp 

m 



2ik 

V 



r(0) 



-1. 



(6.21) 



where fc is the wave vector of the incoming plane wave, 
and ^ is a real parameter that depends on the po- 
tential. The incoming wave is then totally reflected at 



zero energy. There is only one exception to this general 
low-energy behavior, arising for potentials with i/ = 0. In 
that case, the limit fc —s- is different and characterized 
by 



t(0) =to = ±1, r(0) = To 



0, 



(6.22) 



corresponding to a reflectionless situation at zero energy. 
For very small but flnite Eq. (|6.21|) still holds, but 
the reflection and transmission amplitudes change on the 
scale of unity in a narrow momentum region k ^ v. For 
h > V (but still smaller than all other momentum scales) , 
Eq. H6.22|l is reached, and the potential is effectively re- 
flectionless. We now show that for = 0, the conditions 
(i), (ii) and (iii) are satisfled and equivalent to each other, 
while for flnite strictly speaking, none of them holds. 

We first observe from the general definitions of a\ and 
6i, see Eqs. and (jHHII, that Eq. (jlT^ leads to 

the divergence of a\ {b\) for to = 1 (to = —1)- For 7^ 0, 
a\ and b\ are finite so that (i) holds iff (iii) is true. The 
small-j/ regime is also interesting to investigate. Using 
Eq. (Eini), we find that for t{k ~ j/) 1 (-1), ai (61) 
diverges as —2/v. To see if there is a quasi-bound state 
at zero energy, we consider the asymptotic form of the 
wavefunction 



„ikx 



r{k)e 



-ikx 



t{k) e 



ikx 



-00 
-00 



(6.23) 



for which we take fc = 0. If 7^ 0, Eq. (|6.21|l gives V' = 
in the asymptotic region. Since a quasi-bound state at 
zero energy can only be an extended state, we conclude 
that tp vanishes identically in this case. Conversely for 
1/ = 0, we obtain ^p{x) = 1 for x —00 and ip{x) = to = 
±1 for X ^ +00. This means that there is an extended 
quasi-bound state at zero energy. Its parity follows from 
the asymptotic behavior, and we find that it is even if 
to ~ +1 but odd for to = —1. This shows that (ii) and 
(iii), and hence also (i), are equivalent. 

In our specific case, the zero energy mode is found by 
taking u ^ in Eqs. (|6.9|1 and (|6.16|) . This leads to 
(again, e = sgn(a:;)) 



foix) 



_/4(2. 
2e (e 



for bosons, 
for fermions, 



(6.24) 



or Eqs. H6.12|l and (|6.17|) in momentum space. Now fo{x) 
is an even (odd) function for bosons (fermions), corre- 
sponding to a divergence of ai {bi ) , as expected from our 
discussion above. 

To conclude, let us discuss the position of the quasi- 
bound state for large ai or bi by analyzing the poles of the 
scattering matrix S in Eq. I|6.10|l . If a pole appears for the 
scattering amplitude e*''^'^''^ (e*''"*^''^), the corresponding 
bound state is even (odd). For large values of ai, we can 
use the expansion (|6.11|) . and the equation for the even 
bound state is given by 



1 + ikai ~ 0. 



(6.25) 



19 



Analytic continuation to the physical sheet k = iK (k > 
0) gives a real bound state with k = 1/ai for ai > 0, 
with energy = —fi^/mQa\. However, for ai < 0, we 
find, by analytic continuation to the unphysical sheet k = 
virtual bound state with positive energy Ee = 
h^/moal. For the three-boson problem in the BCS limit, 
<iad —f—oo. Consequently, the virtual quasi-bound state 
here comes down to zero from positive energies, and thus 
never becomes a real bound state of the three bosons. 
It simply reaches zero energy in the BCS limit, leading 
to the divergence of ai. The situation is similar for an 
odd bound state. With the expansion Eq. (|6.11|) . we 
find the bound-state equation 1 — ikbi = 0, so that the 
same conclusions as above can be formulated, with just 
a sign difference for bi. Concerning the threc-fermion 
problem, the same scenario for the virtual bound state is 
encountered, since there bad ^ +oo in the BCS limit. 

VII. CONCLUSIONS 



the scattering is described by the three-particle Bethe 
Ansatz equations, i.e., we are dealing with two-body con- 
tact interactions. Then atom-dimer scattering cannot be 
viewed as potential scattering anymore. Although an ef- 
fective potential can be defined, it is non-local and not 
short-ranged . Instead atom-dimer scattering is found 
to be reflectionless: aad diverges for bosons, while bad di- 
verges for fermions. In fact the Bethe Ansatz equations 
remain applicable also for the A^-body problem in the 
BCS limit. 

We have found a novel confined-induced trimcr state 
for bosons. We trace the trimer state numerically from 
the dimer to the BCS limit provided Efimov physics can 
be neglected. Specifically in the dimer limit, we have 
discussed the interplay between this confined trimer state 
and the usual Efimov bound states. The trimer state is 
unique, its energy is nearly universal, and it matches the 
known Bethe Ansatz thrce-particlc bound state in the 
BCS limit. 



In this paper, we presented the results of our study of 
the three-body problem in a quasi-lD confinement. We 
define and calculate two different ID atom-dimer scat- 
tering lengths, ttad and bad, which are directly accessible 
in scattering experiments. Physically, Uad reflects the 
low-energy properties of the scattering phase shift for a 
symmetric wave while bad is the equivalent for an an- 
tisymmetric wave. Both can be inferred from the en- 
ergy dependence of the transmission and the reflection 
coefficients. Technically, we derive a system of integral 
equations for the ID scattering amplitudes with differ- 
ent transverse channel indices. For the sake of simplicity, 
the above system is projected onto the transverse ground 
state. We solve the ground-state equation numerically in 
general and analytically in the dimer and BCS limits. 
Note that the dimer and the BCS limits each correspond 
to different approaches to the 3D problem (|a| ^ a±). 
Therefore we have investigated in detail the role of the 
higher transverse channels in both cases. 

In the dimer limit, indeed it turns out that the higher 
channels contribute on the same footing. Their role is 
described by the 3D equations derived previously |35| . 
Note that the definition of the ID scattering length re- 
mains distinct from that of the 3D scattering length even 
in the dimer limit. We establish a simple analytic rela- 
tion between the latter and aad, valid both for bosons and 
fermions. In this limit bad ^ 0, indicating that atom- 
dimer scattering can be regarded as potential scattering 
with a short-range effective potential. The resulting Uad 
is negative for fermions but positive for bosons pointing 
to the existence of the confinement-induced trimcr state 
in the latter case. Indeed we find such a state from the 
bound-state equation for bosons. Its energy is indepen- 
dent of the bosonic rcgularization parameter R* in the 
large R* limit. 

The BCS limit is even more interesting: we show that 
higher channels can always be neglected here as far as 
the low-energy scattering is concerned. We find that 
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APPENDIX A 

We show here how to obtain the integral representation 
Eq. H1.6|) from the more standard expression (see Ref. 
for example) 

C(l/2,17)= lim (E , J -2iN + n)'/A. 

(AI) 

We define An = (Er^=o(" + ^y^^^) - 2VN + n such 

that limjv_»_(-oo = C(l/2i^^)- Using the integral rep- 
resentations 



1 



dt 



-{n+n)t 



y^n+n Jo V7rt 

dt l-e-(A^+f2)* 



Jq 



t 



N -nt 



and the geometrical summation X]ri=o 
e~(^+-'^'*)/(l — e"*). An can be written as 

An 



^^^Vl-e-* t 



dt /I 



TTt \t 1 - e- 



(1 



(A2) 



-{N+n)t 



The second term in Eq. ljA2p vanishes for large N so that 
we obtain the integral representation of Eq. (|1.6|) . 
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